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ABSTRACT: Long noncoding RNAs (IncRNAs) are perva- 
sively transcribed in the human genome. Recent studies 
suggest that the involvement of IncRNAs in human diseases 
could be far more prevalent than previously appreciated. Here 
we have identified a IncRNA termed Lnc_bc060912 whose 


Cytoplasm 


` 


expression is increased in human lung and other tumors. SS 

Lnc_bc060912 is 1.2 kb in length and is composed of two ( 

exons. The expression of Lnc_bc060912 was repressed by p53. RNA protein* Qf S Nucleus 
Lnc_bc060912 suppressed cell apoptosis. Using a recently a 


developed method for RNA-pulldown with formaldehyde mani apoptosis 

cross-linking, we found that Lnc_bc060912 interacted with 

the two DNA damage repair proteins PARP1 and NPM1. Together, these results suggest that Lnc_bc060912, via PARP1 and 
NPM1, affects cell apoptosis and may play important roles in tumorigenesis and cancer progression. 


WW: the advent and application of high throughput stasis, dosage compensation, and genomic imprinting.’*'* A 


technologies, including microarray and RNA sequenc- number of studies have demonstrated that IncRNAs are 

ing, we have markedly expanded our knowledge of the involved in a wide range of diseases, including cancers,’* 
transcriptome. Transcriptome analyses have revealed that cardiovascular diseases,'* and neurodegeneration diseases.’ 
approximately two-thirds of the mammalian genomic DNA is HOTAIR, one of the most well-characterized IncRNAs, is a 2.2 
transcribed, in contrast to the less than 2% that codes for kb transcript originating from the HoxC locus, and represses 
proteins. ™? Furthermore, organismal complexity correlated the expression of HoxD genes by recruiting the polycomb 
better with the amount of non-protein-coding sequences in the repressive complex 2 (PRC2) to alter the chromatin state of the 
genome rather than with the number of coding genes.** This HoxD locus.’ Moreover, the expression of HOTAIR is 
suggests that noncoding RNAs (ncRNAs) are of significance in increased in primary breast tumors and correlates negatively 
eukaryote cell functions. with the prognosis of the patients." The IncRNA SChLAP1 is 

Besides a number of RNA families with known functions overexpressed in malignant prostate cancers, interacts directly 
(e.g, rRNA, tRNAs, snoRNA, miRNA, etc.), there are tens of with the SWI/SNE complex to remodel chromatin, and serves 
thousands of noncoding transcripts longer than 200 nucleotides as a biomarker for prostate cancer metastasis and mortality.” 
that are commonly referred to as long noncoding RNAs The lincRNA-p21 is regulated by p53, represses cell apoptosis, 
(IncRNAs).° In a comparison with protein-coding genes, and serves as a bridge between p53 and a number of p53 
IncRNAs generally display low abundance and low conservation repressed genes by recruiting hnRNP-K to their genomic loci.'* 
between species and tend to be tissue specific.” Thus, for a LincRNA-p21 was later found to be involved in colorectal 
very long time, IncRNAs were considered to be transcriptional carcinoma."” The LncRNADisease database has thus far 
noise.’ Only in recent decades have IncRNAs received collected and curated more than 1000 IncRNA-disease entries 
considerable attention, as they have been found to participate, and 475 IncRNA interaction entries, including 321 IncRNAs 
through distinct mechanisms at the transcriptional, post- 
transcriptional, or epigenetic level, in a diverse range of Received: March 10, 2015 
biological processes such as stem cell pluripotency, cell Revised: April 6, 2015 
differentiation, cell proliferation and apoptosis, cancer meta- Published: April 7, 2015 
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and 221 diseases from ~500 publications.”° The above studies 
show that IncRNAs may help researchers to understand 
diseases and may have potential clinical significance in disease 
diagnosis, treatment, and prognosis. 

In the present study, we characterized a 1.2 kb long 
noncoding RNA termed Lnc_bc060912. This IncRNA is 
composed of two exons and displayed increased expression in 
tumors compared to adjacent normal tissues, in particular lung 
carcinoma. Moreover, Lnc_bc060912 suppresses cell apoptosis. 
The expression of Lnc_bc060912 is repressed by the expression 
of p53 induced by DNA damage. Finally, we demonstrated that 
Lnc_bc060912 physically interacts with the DNA damage 
repair proteins PARP1 and NPM1. 


M@ MATERIALS AND METHODS 


Cell Lines and Culture Conditions. The H1299, A549, 
and 293T cell lines were used in this study. Among these, 
H1299 and AS49 are lung cancer epithelial cell lines, and 293T 
is an embryonic kidney epithelial cell line. Cell lines were 
cultured using standard protocols provided by ATCC. Briefly, 
cells were cultured at 37 °C in 5% (v/v) CO, incubator. H1299 
and A549 were cultured in RPMI 1640 medium (Life 
Technologies, 11875-093) with 10% fetal bovine serum 
(FBS) (Life Technologies, 16000-044). Cell line 293T was 
cultured in DMEM medium (Life Technologies, DMEM, 
11995-065) with 10% FBS. 

Cells were treated with doxorubicin (dox) in a range of 
concentrations. After 24 h, a Western blot was performed to 
detect the expression level of p53. 

RNA Extraction and Quantitative Reverse Tran- 
scription-Polymerase Chain Reaction (qRT-PCR). Total 
RNA was extracted from cultured cells or nuclear and cytosolic 
fractions with TRIzol regent (Invitrogen, 15596-026) according 
to the manufacturer’s instructions, and then treated with 
recombinant DNase I (Ambion AM2235) enzyme. A 50 ng 
portion of RNA was used for each qRT-PCR by TransScript II 
Green One-Step qRT-PCR SuperMix (TransGen Biotech, 
AQ311-01) according to the manufacturer’s instructions. The 
results were normalized to the expression of GAPDH, and the 
relative expression was calculated by the 6—6 Ct method. The 
qRT-PCR and data collection were carried out on Rotor-Gene 
Q real-time PCR cycler (Qiagen, 9001630). The primer 
sequences are listed in Supporting Information Table 1. 

5’ and 3’ RACE. S’ and 3’ RACE were performed using the 
RLM RACE kit (Ambion AM1700) as recommended by the 
manufacturer. RACE PCR products were obtained using 
GoTaq DNA polymerase (Promega, M3001); the supplied 
GeneRacer primers and the appropriate gene-specific primers 
are summarized in Supporting Information Table 1. RACE 
PCR products were separated on a 1.5% agarose gel. Gel 
products were extracted with the Gel and PCR Clean-Up 
System (Promega, A9282), cloned into the pGEM-T Vector 
Systems I (Progema, A3600), and sequenced bidirectionally 
using the M13 forward and reverse primers by Sanger 
sequencing at Invitrogen. At least five colonies were sequenced 
for every RACE PCR product that was gel purified. 

Subcellular Fractionation. The separation of the nuclear 
and cytosolic fractions was performed using the NE-PER 
Nuclear and Cytoplasmic Extraction Reagents (Thermo 
Scientific, 78833) according to the manufacturer’s instructions. 
RNA was extracted from the two components. Then, qRT-PCR 
was performed to assess the relative proportion of specific RNA 
molecules in the nuclear and cytoplasmic fractions. GAPDH 
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and f-actin serve as the cytosolic controls, and U1 snRNA 
serves as the nuclear control. The primer sequences are listed in 
Supporting Information Table 1. 

Northern Blot. Northern blots were carried out using the 
DIG Northern Starter Kit (Roche, 12039672910) according to 
the manufacturer’s instructions. A 30 ug portion of total RNA 
and a 3 pg portion of total RNA were loaded on 0.8% 
denaturating formaldehyde agarose gels for detecting IncRNA 
Lnc_bc060912 and mRNA GAPDH, respectively. Digoxigenin 
(DIG)-labeled antisense probes were synthesized by in vitro 
transcription using T7 RNA Polymerase (Roche, 
10881767001) and DIG RNA Labeling Mix (Roche, 
11277073910) and taking a purified DNA PCR product with 
T7 promoter as template. Sequences of the PCR primers for 
probe templates are listed in Supporting Information Table 1. 

siRNA Mediated Knockdown. siRNAs used to knock 
down target IncRNAs or mRNAs were designed by Custom 
RNAi Design Tool on the Integrated DNA Technologies 
(IDT) website and synthesized by the GenePharma company. 
A nonspecific siRNA used as negative control (NC) was also 
purchased from the GenePharma company. siRNA sequences 
(S’ to 3’) for the knockdown experiments are listed in 
Supporting Information Table 1. Cells were seeded in 6-well 
plates at a desired concentration and transfected with 100 pmol 
of siRNA oligonucleotide after 8 h plating. The transfection was 
performed with Lipofectamine 2000 (Invitrogen, 11668-019) 
according to the manufacturer’s instructions. 

Cloning. All PCR primers used for cloning are listed in 
Supporting Information Table 1. The DNA polymerase for 
cloning was PfuUltra II Fusion HS DNA Polymerase (Agligent, 
600670), and restriction endonucleases were from New 
England Biolabs (NEB). The full-length Lnc_bc060912 
transcript and TPS3 open reading frame were amplified from 
AS49 cDNA and cloned into pcDNA3.0 plasmid and 
p3xFLAG-CMV-10 plasmid. Anti-pS3 (Santa Cruz, sc-126) 
and antiflag (Sigma, F3165) were used. 

Overexpression. Overexpression of ectopic transcript 
sequences was performed by transient transfection using a 
plasmid containing the transcript sequence and blank plasmid 
as negative control. Transfection reagent was the same as that 
used in siRNA to mediate knockdown experiments. 

Cell Proliferation Assay. Cell proliferation was monitored 
using the CellTiter 96 Aqueous One Solution Cell Proliferation 
Assay (MTS) (Promega, G3582). The transfected cells were 
plated in 96-well plates (2000 cells/well). Cell proliferation was 
determined every 24 h according to the manufacturer’s 
instructions. 

5-Ethynyl-2’-deoxyuridine (EdU) Incorporation Assay. 
Upon 48 h after transfection, the 5-ethynyl-2’-deoxyuridine 
(EdU) incorporation assay was performed with the Click-iT 
EdU Alexa Fluor 488 Imaging Kit from Life Technologies 
(C10337) according to the manufacturer’s instructions. 

Cell Apoptosis Assay and Cleaved Caspase-3 Meas- 
urement. For cell apoptosis quantification, cell viability was 
determined by staining with Annexin V, Alexa Fluor 488 
conjugate (Life Technologies, A13201), and propidium iodide 
(Life Technologies, P21493) and detected by flow cytometry 
(BD FACSVERSE) according to the Annexin V flow cytometry 
protocol. Cleaved caspase-3 was measured with Western blot 
using the Cleaved Caspase-3 antibody (ABGENT, AP3725a). 
The cells were collected for flow cytometry assays or total 
protein lysate 48 h after transfection. 
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Figure 1. Identification and characterization of the Lnc_bc060912 as a long noncoding RNA. (A) Transcription profiles of Inc_bc060912 in four 
types of paired tumor tissue samples by microarray. N and C represent the normal and cancer tissues, respectively. (B) Genomic location of 
Lnc_bc060912 relative to upstream and downstream genes. (C) Northern blots of the Lnc_bc060912 and GAPDH transcripts. (D) Quantification 
of Lnc_bc060912 expression by qRT-PCR of fractionated A549 cell lysates. U1 RNA serves as a positive control for nuclear gene expression. 
GAPDH and /-actin serve as controls for cytoplasmic gene expression. Bars indicate SEM, n = 3. (E) Relative expression of Lnc_bc060912 and its 
neighbor genes after knockdown of Lnc_bc060912 by four different siRNAs. Bars indicate SEM, n = 3, *P < 0.05, **P < 0.01. 


RNA-Pulldown. RNA-pulldown was carried out as 
previously described’! with some modifications. Antisense 
DNA probes were designed against a Lnc_bc060912 full-length 
sequence using the online design software at the Stellaris Fish 
Probe website (http:/ /www.singlemoleculefish.com). A set of 
probes against LacZ RNA was also generated as controls. There 
were 10 probes generated against Lnc_bc060912 and 8 for 
LacZ. All probes were biotinylated at the 3’ end with an 18- 
carbon spacer arm (Invitrogen). Originally 10 plates of 15 cm 
dish 293T cells were cross-linked by 1% formaldehyde. Cross- 
linked cells were lysed in lysis buffer at 100 mg/mL (SO mM 
Tris pH 7.0, 10 mM EDTA, 1% SDS, add DTT, PMSF, 
protease inhibitor, and RNase inhibitor before use) on ice for 
10 min, and sonicated using Ultrasonic Cell Crusher on ice 
until the lysate became transparent and clear. The cell lysate 
was separated into two equal aliquots, one for probes targeting 
the IncRNA, and the other for probes targeting LacZ as 
controls. For every aliquot, we used 600 pmol probes and 200 
HL streptavidin magnetic Cl beads. Probes and beads were 
incubated with cell lysate, and the beads were washed 5 times 
with washing buffer (2xSSC, 0.5% SDS, add DTT and PMSF 
fresh). We separated beads into two parts, 1/20 for RNA 
elution and 19/20 for protein elution. For elution of the RNA, 
beads were treated with Proteinase K, and RNA was extracted 
with Trizol reagent. Then, qRT-PCR was performed to 
determine the enrichment of target RNA compared to that of 
the control. For elution of the protein, cross-linking was 
reversed by incubating beads with -mercaptoethanol at a final 
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concentration of 2.5% in 1x NuPAGE LDS Sample Buffer (Life 
Technologies, NP0007) at 96 °C for 30 min. The proteins were 
separated in a 4—12% Bis-Tris Gel (Life Technologies, 
NP033SBOX), stained with the Silver Staining Kit (Life 
Technologies, LC6070), and followed by mass spectrometry 
(MS) identification. 

RNA Immunoprecipitation (RIP). The RIP experiments 
were performed using a Magna RIP RNA-Binding Protein 
Immunoprecipitation Kit (Millipore, 17-700) according to the 
manufacturer’s instructions. The RIP product (1/10) and 1% of 
the cell lysate input were used for Western blot as described 
above, while the rest of the RIP product was used for RIP qRT- 
PCR. The antibodies for RIP assays were as follows: normal 
rabbit IgG (Santa Cruz, sc-2027), normal mouse IgG (Santa 
Cruz, sc-2025), anti-PARP-1 antibody (Santa Cruz, sc-8007), 
and anti-nucleophosmin antibody (Abcam, ab52644). All RIP 
assays were performed in biological duplicate. 


E RESULTS 


Identification and Characterization of the Long 
Noncoding RNA Lnc_bc060912. We collected 76 x 2 
paired tumor and adjacent normal tissue samples, including 20 
X 2 paired gaster samples, 20 X 2 paired colon samples, 16 X 2 
paired liver samples, and 20 X 2 paired lung samples. For the 
identification of IncRNAs that were differently expressed in 
tumors, total RNA was extracted from the samples and 
hybridized to the mRNA + IncRNA V2.0 microarray on an 
Agilent platform. Through a pipeline for analysis of the 
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expression difference, we identified a number of significantly 
differentially expressed IncRNAs from all the tissue types 
examined (unpublished data). There is one transcript, 
BC069212, that showed increased expression in several tumors 
compared to paired normal tissues, and displayed the most 
significant difference in lung cancer (Figure 1A). Subsequently, 
we attempted to characterize BC069212 in lung cancer cell 
lines. 

BC069212 is located in the chr22 q13.31 region between 
gene GTSE1 and gene TRMU (Figure 1B). It has been 
annotated as BC069212 in UCSC Genes. To determine both 
ends of this transcript, we performed S’ and 3’ RACE assays in 
AS49 cells. We discovered that the transcript is polyadenylated, 
contains two exons, and maps to chr22:46727009—46729596, 
and the intron is located at chr22:46727630—46728978. 
Furthermore, the Northern blot was also performed to confirm 
the 1.2 kb full length of the transcript (Figure 1C). Next, we 
used the CPAT algorithm” to assess the coding potential of 
the transcript. We obtained a coding probability of 0.1076 
which strongly suggests that it is a noncoding RNA. Then, 
isolation of the nuclear and cytosolic fractions of A549 cells 
showed that the transcript was mainly located in the nucleus 
(similar to that of Ul, Figure 1D). Taken together, the 
transcript is a nuclear-retained long noncoding RNA, and we 
termed the transcript as Lnc_bc069212. 

As there are a number of IncRNAs reported to affect the 
neighboring protein-coding gene expression, ^ we examined 
the gene expression of GTSEl and TRMU to determine 
whether Lnc_bc060912 exerts its function in cis. We 
consequently knocked down the expression of Lnc_bc060912 
by RNA interference using four different siRNAs, and examined 
the mRNA expression level of GTSE1 and TRMU. Neither of 
the genes was significantly affected by the knockdown of 
Lnc_bc060912 (Figure 1E), and it is thus not likely that 
Lnc_bc060912 acts in cis to regulate the expression of its 
nearest neighbors. 

Lnc_bc060912 Represses Cell Apoptosis. Since 
Lnc_bc060912 showed higher expression in tumor samples 
than adjacent normal tissues did, we next aim to explore its 
function during the tumorigenesis. From above, we selected 
two more effective siRNAs to knockdown Lnc_bc060912 in 
AS49 cells. MTS assay showed that knockdown of 
Lnc_bc060912 significantly impaired cell viability in the A549 
cell line (Figure 2A). Further, we assayed the incorporation of 
5-ethynyl-2’-deoxyuridine (EdU) of DNA in cells by 
fluorescence-activated cell sorting (FACS). Also, we observed 
reduced EdU incorporation in cells treated with siRNAs 
compared to controls (Figure 2B,C). These results suggested 
that Lnc_bc060912 plays a role in regulating cell proliferation. 

A reduction in cell proliferation will likely be a result of an 
increase in cell apoptosis. We then examined the effect of 
Lnc_bc060912 knockdown on apoptosis. To this end, we 
assayed the proportion of the cell population undergoing 
apoptosis by measuring Annexin-V by FACS analysis (Figure 
2D). We observed a significant increase in the number of 
apoptotic cells after knockdown of Lnc_bc060912 relative to 
the siRNA control (Figure 2E). We also observed an increase in 
caspase-3 cleavage after knockdown of Lnc_bc060912 relative 
to controls (Figure 2F). Thereafter, we performed an 
Lnc_bc060912 overexpression assay in the A549 cell line by 
transfecting a plasmid containing the Lnc_bc060912 sequence. 
The expression level of Lnc_bc060912 in the transfected cells 
was over 400 times greater than that of control cells 
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Figure 2. Lnc_bc060912 represses cell apoptosis. (A) The absorbance 
on OD 490 nm in MTS assay every 24 h after transfection of siRNAs. 
(B, C) EdU incorporation ratio after Lnc_bc060912 knockdown. (B) 
Readout from the flow cytometer. The M2 area indicates cells 
incorporating EdU (i.e., synthesizing DNA). The percentage of cells in 
each area is indicated. (C) Quantification of results shown in part B. 
EdU incorporation ratio of three independent replicate experiments. 
Bars indicate SEM, n = 3, **P < 0.01. (D, E) Cell apoptosis ratio after 
Lnc_bc060912 knockdown. (D) Readout from the flow cytometer. 
The x axis represents Annexin-V, and the y axis represents PI staining. 
LR indicates early apoptotic cells. UR indicates terminal apoptotic 
cells. The percentage of cells in LR and UR is indicated. (E) 
Quantification of the results in part D. Percentage of apoptotic cells, 
with average of three independent replicate experiments. Bars indicate 
SEM, n = 3, *P < 0.05, **P < 0.01. (F) Increase in caspase-3 cleavage 
48 h after Lnc_bc060912 knockdown by siRNA transfection. (G, H) 
Cell apoptosis after Lnc_bc060912 overexpression. (G) Read out from 
the flow cytometer. (H) Quantification of the result in part G. 
Percentage apoptotic cells, average of three independent replicate 
experiments. 


(Supporting Information Figure S1), and we observed a low 
apoptosis rate in the cells overexpressing Lnc_bc060912 
(Figure 2G,H). 
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Figure 3. p53 suppresses Lnc_bc060912 expression. (A) Effect of doxorubicin (dox) on p53 in A549 cells (Western blot). (B) Effect of dox on 
Lnc_bc060912 and GTSE1 expression in A549 and H1299 cells. (C) Demonstration of ectopic p53 overexpression level in H1299 cells (Western 
blot). (D) Effect of p53 overexpression on Lnc_bc060912 and GTSE1 in H1299 cells. (E) p53 knockdown by siRNAs in A549 cells (Western blot). 
(F) Effect of p53 knockdown on Lnc_bc060912 expression in A549 cells. Bars indicate SEM, n = 3, *P < 0.05, **P < 0.01, ***P < 0.001. 


p53 Suppresses the Expression of Lnc_bc060912. The 
tumor suppressor p53 is known to be involved in apoptosis, 
and we therefore investigated the relationship between 
Lnc_bc060912 and p53. We treated A549 cells (expressing 
wild type p53) and H1299 cells (in which the p53 gene has 
been homozygously deleted) with doxorubicin (dox), and 
found that p53 was upregulated upon dox treatment in a dose- 
dependent manner as expected in the A549 cell line (Figure 
3A). We also performed qRT-PCR to examine the transcrip- 
tional levels of Lnc_bc060912 and its upstream gene GTSE1. 
The results showed that, in the A549 cell line, the RNA levels 
of Lnc_bc060912 and GTSE1 decreased remarkably with the 
increase in dox concentration, while the expression levels of 
neither gene displayed considerable decrease in the H1299 cell 
line compared with that in the AS49 cell line under the same 
treatment (Figure 3B). 

We wondered whether increasing the level of p53 in H1299 
would suppress the expression of Lnc_bc060912 and GTSE1 in 
these cells. To this end we constructed a p53-p3xflag fusion 
plasmid, and transfected it into H1299 cells. After 48 h, we 
examined the expression of p53 by Western blot and the 
expression of Lnc_bc060912 and GTSE1 by qRT-PCR (Figure 
3C,D). The enhanced expression of p53 markedly decreased 
the expression of Lnc_bc060912 and its nearby gene, GTSE1 
(Figure 3D). Previous microarray analysis showed that the 
expression of Lnc_bc060912 and GTSE1 were strongly 
correlated (Supporting Information Figure $2). These implied 
that Lnc_bc060912 and GTSE1 were influenced by a common 
regulator (i.e. p53, or one of its downstream genes). 

Thereafter, we knocked down p53 in the A549 cell line using 
three different siRNAs, which all effectively reduced the 
expression of p53 (Figure 3E,F). The expression level of 
Lnc_bc060912 showed a significant increase in the presence of 
a lower p53 level (Figure 3F). In conclusion, the tumor 
suppressor p53 represses the expression of Lnc_bc060912. 


2899 


Lnc_bc060912 Physically Interacts with PARP1 and 
NPM1. Having characterized the function of Lnc_bc060912, 
we also wanted to explore in further detail the mechanism by 
which Lnc_bc060912 mediates apoptosis. RNA-pulldown 
coupled with mass spectrometry has provided a method to 
identify the potential interactions between RNAs and 
proteins. Using the traditional RNA-pulldown method 
which employed sense and antisense Lnc_bc060912 tran- 
scripts, we observed no different protein bands (Supporting 
Information Figure S3). We then try to use a ChIRP 
experiment coupled with mass spectrometry to discover the 
proteins which specifically interact with Lnc_bc060912. 

By using a ChIRP experiment in the 293T cell line, we 
enriched endogenous Lnc_bc060912 above 1000 times by 
Lnc_bc060912 antisense probes compared to LacZ antisense 
probes (Figure 4A), observed two specific protein bands, and 
subjected them to mass spectrometry (Figure 4B). We obtained 
the two proteins poly(ADP-ribose) polymerase 1 (PARP1) and 
nucleophosmin (NPM1), from each of the two respective 
bands. 

To further confirm the interactions between Lnc_bc060912 
and PARP1, and Lnc_bc060912 and NPM1, we performed 
RNA immunoprecipitation (RIP) with antibodies against 
PARP1 and NPM1. Both antibodies worked well, and the 
precipitates were significantly enriched for Lnc_bc060912 
(Figure 4C,D). Thus, the results suggested the existence of 
physical interactions between Lnc_bc060912 and PARP1, and 
between Lnc_bc060912 and NPM1. 

PARP! is located in the nucleus and localizes at sites of DNA 
damage. NPM1 is located in the nucleolus, but it can be 
translocated to the nucleoplasm in some cases. The nuclear 
localization of the two proteins is in accordance with that of 
Lnc_bc060912 (Figure 1E). PARP1 and NPM1 are both 
involved in DNA damage repair, °° and have an effect in 
preventing apoptosis, which is also consistent with the observed 
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Figure 4. Lnc_bc060912 interacts with PARP1 and NPMI1. (A) 
Pulldown assay enriches Lnc_bc060912 RNA more than 1000-fold. All 
relative abundances were compared to 1% input. (B) SDS-PAGE gel of 
proteins pulled down by the Lnc_bc060912 and LacZ probes. Bands 
indicated by arrows were submitted for mass spectrometry. (C) RNA 
immunoprecipitation (RIP) of PARP1 enriches Lnc_bc060912. U1 
and GAPDH are negative controls. All relative abundances were 
compared to 1% input. (D) RNA immunoprecipitation (RIP) of 
NPMI1 enriches Lnc_bc060912. Ul and GAPDH are negative 
controls. All relative abundances were compared to 1% input. All 
experiments were performed in three independent replicates. Bars 
indicate SEM. 


activity of Lnc_bc060912. Consequently, the increased 
expression of Lnc_bc060912 in tumors is in accordance with 
the role in preventing cell apoptosis through interactions with 
PARP1 and NPM1. 


M DISCUSSION 


It has become increasingly clear that mammalian genomes 
transcribe thousands of IncRNAs.”’ °° Dysregulation of 
IncRNAs is often associated with human diseases, in particular 
cancer,°' and recent progress suggests that the involvement of 
IncRNAs in human diseases could be far more prevalent than 
previously appreciated.* 

In the current study, we have identified a novel endogenous 
IncRNA Lnc_bc060912 which is highly expressed in four tumor 
types (gaster, colon, liver, and lung) compared with adjacent 
normal tissues. Human Body Map Project data shows that 
Lnc_bc060912 is expressed in nearly all tissues, with slight 
differences in expression (Supporting Information Figure 4)? 
Thus, expression of the Lnc_bc060912 RNA may be required 
to maintain normal cell functions in the majority of the body 
tissues. 

LncRNAs participate in a wide repertoire of biological 
processes, and almost every step in the life cycle of genes, from 
transcription to post-transcription, can be influenced by 
IncRNAs.””!'3! Here, we have shown that IncRNA 
Lnc_bc060912 acts to repress cell apoptosis. This function is 
in accordance with the oncogenic function of Lnc_bc060912, 
indicated by its increased expression in tumors. A number of 
studies have reported that dysregulation of IncRNAs affects cell 
function. For example, the ncRNA SChLAP1, which is highly 
expressed in malignant prostate tumors, promotes cancer cell 
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invasiveness and metastasis.” The IncRNA PINT is down- 
regulated in primary colorectal tumors, and its overexpression 
inhibits the proliferation of tumor cells.” 

The fact that the transcription factor p53 acts upstream of 
Lnc_bc060912 in the regulatory network, and negatively 
regulates the expression of Lnc_bc060912, is also an interesting 
observation. TP53 plays an important role in apoptosis, 
genome stability, and inhibition of angiogenesis in multicellular 
organisms.** Huarte and colleagues reported that numerous 
lincRNAs are activated in a pS3-dependent manner, and 
characterized one lincRNA, lincRNA-p21, which serves as a 
transcriptional repressor in the p53 pathway and plays a role in 
triggering apoptosis." Zhang and colleagues reported that the 
expression of the p53-induced IncRNA TUGI serves as an 
independent predictor for overall cancer survival, and that 
TUG1 knockdown significantly promoted cell proliferation.*° 
In contrast, Lnc_bc060912 is suppressed by p53, which 
suggests that the IncRNAs suppressed by p53 may also be 
important for cell function and tumorigenesis. Also, unlike 
lincRNA-p21, Lnc_bc060912 represses cell apoptosis. 

LncRNAs affect cellular functions through various mecha- 
nisms, such as interactions with chromatin modification 
proteins, chromatin remodeling proteins, transcription factors, 
and microRNAs. The polycomb repressive complex 2 (PRC2) 
catalyzes histone H3 on lysine 27 trimethylation for epigenetic 
silencing. Researchers have discovered that a number of 
IncRNAs, such as HOTAIR,’® PINT,” and TUG1,” interacted 
directly with PRC2 and guided PRC2 to target genes for 
H3K27 trimethylation. Chromatin remodeling complexes such 
as SWI/SNEF complex and transcription factors such as c-Myc 
have been observed to interact with IncRNAs.'”*° In the 
present study, we observed that Lnc_bc060912 interacted 
physically with PARP1 and NPM1, neither of which having 
previously been reported to interact with IncRNAs. PARP] is 
involved in the regulation of various important cellular 
processes such as differentiation, proliferation, and tumor 
transformation and also in the regulation of the molecular 
events involved in the recovery of cell from DNA damage.” 
NPM1I has multiple functions in genome stability, DNA repair, 
and regulation of apoptosis.”° The functions of the two proteins 
are in accordance with the observations made here concerning 
the functions of Lnc_bc060912; thus, it makes sense to assume 
that Lnc_bc060912 represses apoptosis through its interaction 
with PARP1 and NPM1. One report suggests that PARP1 and 
NPM interact with each other and accumulate in transcrip- 
tionally active nucleoli;*” however, this question of these two 
proteins interacting together or independently with 
Lnc_bc060912 to repress apoptosis still needs to be explored. 

In addition to elucidating some functional aspects of 
Lnc_bc060912 IncRNA, our study also represents a certain 
methodological improvement over current methods. Currently 
the most commonly used method to explore proteins 
interacting with IncRNAs is the RNA-pulldown procedure 
originally described by John L. Rinn.’* In brief, in this method 
they incubate a full-length biotin-labeled IncRNA and its 
antisense as control with cell lysate, and then subject the 
capture proteins to mass spectrum analysis. The captured 
proteins depend on IncRNA-protein affinity in vitro, which 
results in the binding of many nonspecific proteins. In the 
present study, we used ChIRP to capture proteins interacting 
with the IncRNA.*"* In this method, IncRNAs and proteins 
are cross-linked by formaldehyde to reduce nonspecific 
bindings, and antisense probes were used to capture in vivo 
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IncRNAs. We observed a 1000-fold enrichment in the activity 
of IncRNA Lnc_bc060912 using targeting probes when 
compared to the use of LacZ probes, and identified two 
binding proteins. Subsequently, we used RIP-qPCR to confirm 
the interaction between IncRNA and the two proteins. The 
enrichment of IncRNA Lnc_bc060912 by immunoprecipitation 
of the two proteins was also more than 1000-fold above input. 
Thus, we have provided a new and better method for RNA- 
pulldown. 
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